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We construct minimal electronic models for a newly discovered superconductor LaOi_a;Fi:BiS2 
(Tc = 10. 6K) possessing BiS2 layers based on first principles band calculation. First, we obtain 
a model consisting of two Bi 6p and two S 3p orbitals, which give nearly electron-hole symmetric 
bands. Further focusing on the bands that intersect the Fermi level, we obtain a model with 
two p orbitals. The two bands (per BiS2 layer) have quasi-one-dimensional character with a double 
minimum dispersion, which gives good nesting of the Fermi surface. At around x ~ 0.5 the topology 
of the Fermi surface changes, so that the density of states at the Fermi level becomes large. Possible 
pairing states are discussed. 



For the past several decades, superconductors with 
layered lattice structures such as the cupratesjl|, or- 
ganic conductors pi], MgEzQ, Sr2Ru04[l, Na2,Co02[3, 
(Hf,Zr)NCl[6,], and iron pnictidesQ have attracted much 
attention. These layered superconductors have been of 
interest in many aspects like high and/or unconven- 
tional pairing mechanisms. Quite recently, Mizuguchi et 
al. have discovered superconducting materials that pos- 
sess BiS2 layers, where the Bi and S atoms are aligned 
alternatively on a square lattice. The materials found so 
far are 3140483 [^l and LaOi-,;F:,BiS2'9'l which have Tc 
of 8.6K and 10. 6K, respectively. Fig[l] shows the lattice 
structure of LaOBiS2, where partial replacement of O by 
F (x = 0.5) provides electron doping and gives rise to 
superconductivity. These findings strongly suggest that 
materials with BiS2 layers provide yet another family for 
layered superconductors exhibiting double-digit Tc, and 
it is of special importance to understand the underlying 
electronic structure of these materials. 

Recent theoretical studies on superconductivity have 
shown that effective models constructed from first prin- 
ciples band calculation can provide solid basis for the 
study on the mechanism of superconductivity, especially 
for materials with complicated band structures such as 
the iron pnictides[lo|. In the present paper, we perform 
first principles band calculation of LaOBiS2, from which 
we construct maximally localized Wannier orbitals to ob- 
tain the effective tightbinding models, i.e., the electron 
kinetic energy part of the effective Hamiltonian. The 
models consist of two or four bands that have quasi-one- 
dimensional character, which are hybridized to give a two 
dimensional Fermi surface. At around x ~ 0.5 the topol- 
ogy of the Fermi surface changes, so that the density of 
states at the Fermi level becomes large. We show that the 
quasi-one-dimensional nature results in a nesting of the 
Fermi surface, which gives rise to enhanced irreducible 
susceptibility along the diagonals that intersects the wave 
vectors (0,0) and (tt, tt). We discuss possible supercon- 
ducting states. 

First principles band calculation of the mother 
compound LaOBiS2 is performed using the Wien2K 
package Uj and adopting the lattice structure given in 



ref.jl^]. Space group is P4/nmm and a = 4. 05 A, 
I3.74A. La, S and Bi take the 2b positions with zl 




FIG. 1. The lattice structure of LaOBiS2. 
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0.0936, = 0.813, zs, = 0.379 and zb^ = 0.6263, and 
O takes the 2a position. Here we present results with- 
out the spin-orbit couphng, although this coupling does 
have some effect on the band structure. The calculation 
result is shown in Fig 15] The conduction bands around 
the Fermi level consists mainly of in-plane Bi Qp orbital 
character, mixed with in-plane S 3p. 

From this band calculation, we obtain maximally lo- 
calized Wannier orbitals[l3,[l3]j which enables us to con- 
struct tightbinding models that correctly reproduces the 
original first principles band structure around the Fermi 
level. First, we construct a 24 orbital model, which con- 
sists of six Bi 6p, twelve S 3p and six O 2p orbitals. The 
band structure of the obtained tight binding model is 
shown in FiglSlJa). Here, the thickness of the lines rep- 
resent the weight of the in-plane p orbitals within the 
BiS2 layers. Since the BiS2 layers are without doubt the 
origin of the superconductivity, we can further obtain a 
model that omits the O 2p orbitals, the ip orbitals of the 
out-of-plane S , and also the pz orbital of the in-plane S 
(whose bands lie away from the Fermi level), and we are 
left with an eight orbital model (band not shown). There 
are eight orbitals because there are two BiS2 layers per 
unit cell, and each Bi and in-plane S has two p orbitals. 



By further neglecting the small interlayer coupling be- 
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FIG. 3. (a) The 24 orbital model and (b) the four orbital 
model. In (b), the dashed (dotted-dashed) lines denote the 
Fermi energy for the doping ratio 5 = 0.5 {5 = 0.25). 
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FIG. 4. (a) The tightbinding model on the Bi square lattice 
is shown along with the hoppings of px and py orbitals. The 
stars denote the position of the S atom, whose p orbitals are 
explicitly considered in the four orbital model, but are inte- 
grated out in the two orbital model, (b) The band structure of 
the two orbital model. The dashed (dotted-dashed) lines de- 
note the Fermi energy for the doping ratio 5 — 0.5 {5 = 0.25). 
(c) The Fermi surface of the two orbital model for the doping 
rate of 5 — 0.25 and S = 0.5. (d) A schematic figure of the 
Px band along (— tt, — tt) — ^ (0,0) — >■ (-|-7r,+7r). Three types 
of nesting vectors are shown by the arrows. 
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tween the neighboring BiS2 layers, we end up with a 
two-dimensional "p-p" four orbital model consisting of 
two Bi 6p and two S 3p orbitals. The band structure 
of this model is shown in FiglSl^b). We write the model 
Hamiltonian in the form 

ij a 

+ t{xj - Xi, Vj - y,; V, m)c]^^c,^<t] + ^ e^^w (1) 

where cj^^ creates an electron with spin a on the /i-th or- 
bital in the i-th unit cell, and rii^g- — cl^^Ci^a- Then, the 
parameters of this four orbital model is given in tabl^ 
The two p orbitals are denoted as px and py, where the 
X-Y axis are rotated by 45 degrees from the x-y axes (an 
figure of the four orbital model is not presented, but see 
FiglU^a) for the X and Y axes.) It is interesting to note 
that the upper and the lower bands are roughly symmet- 
ric with respect to the gap. 

Even more simple model is the one which focuses only 
on the bands that intersect the Fermi level. By extract- 
ing these bands using the maximally localized Wannicr 
orbitals centered at the Bi sites and neglecting the intcr- 
layer hoppings, the Hamiltonian reduces to a two dimen- 
sional two orbital model (FigUJ^a)), whose band structure 
is shown in FigUJ^b) and the hopping parameters given 
in tabljin Assuming a rigid band, the Fermi surface for 
this model is obtained in FiglUJc) for two doping ratios 
6 = 0.25 and S = 0.5. Here 5 is defined as the number 
of doped electrons per Bi site, and satisfies 5 = x in an 
ideal situation. At around S ^ 0.5 the topology of the 




FIG. 5. The largest eigenvalue of the irreducible susceptibility 
matrix for S = 0.25 and 0.5. 
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FIG. 6. (a) sign conserving s-wave (b) sign reversing s-wave, 
and (c) d-wave superconducting gaps. 



TABLE I. Hopping parameters t(Ax, Ay; /i, f) or the on- 
site energies for the four orbital model. The orbitals 
= 1,2,3,4 correspond to Bi py, px, S py, px orbitals, 
respectively. We follow the notations in ref . [lOl] . i.e., /, and 
aa corresponds to t{—Ax, — Ay; jjb, v), t{Ay, Ax; fj,, v), respec- 
tively, where '±' means that the corresponding hopping is 
equal to ±t[Ax, Ay; fi, v), respectively. 
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TABLE II. Hopping parameters and on-site energies for the 
two orbital model, = 1, 2 correspond to py, Px, respec- 
tively Here t' = 0.880, t = -0.167, and txy = 0.107 
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Fermi surface changes, so that the density of states at 
the Fermi level becomes large around this band filling. 
The Px and py bands have essentially one-dimensional 
character, where the main hopping integrals (t') exist be- 
tween next nearest neighbor sites, as shown in Fig|4l^a) 
and tableHH Thus the one-dimensional bands essentially 
have the forms 

Sxik) = 2t' cos{kx + ky ), ey{k) — 2t' cos{kx — ky ) (2) 

which has a double well band dispersion along {kx,ky) =^ 
(-TT, -tt) (0,0) (-|-7r,-f7r) or (-7r,+7r) (0,0) 
(+7r, — tt), as shown in Fig|4l^d). The intra-orbital (t) as 
well as the inter-orbital {txy) nearest neighbor hoppings 
give the two dimensionality of the system. 

The one-dimensional nature of the bands provides 
good nesting of the Fermi surface. To see this effect, we 
calculate for the two orbital model the 4x4 irreducible 
susceptibility matrix as k summation of the product of 
the 2x2 Green's function matrix. In FiglSJ we show the 
largest eigenvalue of the irreducible susceptibility matrix 
for the doping ratio of i5 = 0.25 and 5 = 0.5. The diago- 
nal structures that go through (0,0) or (tt, tt) are due to 
the nesting shown by the arrows in Fig|4l[d) . Note that if 
the nearest neighbor hoppings do not exist and the bands 
have the form in eg.©, the nestings (i)~ (iii) are equiv- 
alent, while this equivalence is lost in the presence of t 
(and additional hoppings). 
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Finally, let us discuss the possible pairing states and 
the superconducting gap structure. Since the many-body 
term of the effective Hamiltonian is not determined here, 
there are several possibilities at present. The most rel- 
evant bands have mainly 6p character, which gives wide 
spread of the Wannier orbitals, so that the electron- 
electron interactions may not be very strong and short- 
ranged as in the Sd-orbital based materials such as the 
cuprates and iron-based superconductors. In that case, 
the electron-phonon interaction can be playing the main 
role in the Cooper pairing, and the good nesting of the 
Fermi surface may cooperate to give enhanced attrac- 
tive pairing interaction around the nesting vectors. This, 
combined with the large density of states for x ~ 0.5 
due to the change in the topology of the Fermi surface, 
can give rise to a s-wave pairing with constant gap sign 
as shown in Fig|6l^a). A more exotic possibility for the 
same type of gap is related to the fact that the system is 
close to a band insulator, where electron-hole excitations 
might be playing an important role in the pairing |2 51]. In 
fact, if we consider electron-electron interactions in the 
four orbital model, the system can be viewed as two in- 
teracting charge-transfer-type insulating bands. This is 
similar to the situation studied in ref.[26|, where the oc- 
currence of superconductivity was proposed for systems 
with a metallic band interacting with a charge-transfer- 
type insulating band. In this context, it should be noted 
that the analysis for the Bi 6p-S 3p four orbital model 
can give different results from those for the two orbital 
model regarding the superconducting state, because the 
band filling is different, i.e., the four orbital model is a 
nearly half-filled, electron-hole symmetric system , while 
the two orbital model has small band filling. 

On the other hand, if we assume a short ranged re- 
pulsive interaction, the present two orbital model pro- 
vides a fundamental problem of the pairing state in repul- 
sively interacting quasi-one-dimensional systems, apart 
from its relevance to the BiS2 layers. Namely, the oc- 
currence of superconductivity in repulsively interacting 
one-dimensional systems with a double well type band 
structure was studied in the late 1990's, where the in- 
teraction (ii) in FigHlJd) dominates for certain parame- 
ter regime to give a gap that reverses its sign between 
the inner and the outer Fermi points 17- 23|. The exis- 
tence of t in FigHJa) is essential in this case, because 
this makes the interactions (i)/(iii) and (ii) inequivalent. 
In fact, by adding on-site intra- (U) and inter-orbital in- 
teractions {U' = 2U/3 and the Hund's coupling and the 
pair hoppings J = J' = U/Q), and applying multior- 



bital random phase approximation [15|, |16j to obtain the 



driven by the nestings given in Fig|4l^d). The gap (b) 
indeed exhibits sign reverse between the inner and the 
outer Fermi surfaces. On the other hand, gap (c) is ob- 
tained by reversing the gap sign for both interactions (i) 
and (ii). To make the gap have even parity requires ad- 
ditional nodes at kx = iky which do not intersect the 
Fermi surface for small doping ratio. 

Finally, yet another possibility related with the quasi- 
one-dimensionality is the spin-triplet pairing. In fact, the 
present band structure also has similarity with the quasi- 
one-dimensional bands (the hybridized dxz/yz bands) in 
Sr2Ru04, where some theoretical studies have suggested 
the possibility of spin-triplet pairing originating from 
these bands 2ll424| . Nuclear magnetic resonance experi- 
ments will provide a test for this possibility. 

To summarize, we have obtained effective tightbinding 
models for the newly discovered superconducting BiS2 
layers by performing first principles band calculation for 
LaOBiS2 and exploiting the maximally localized Wan- 
nier orbitals. The model consists of px and py orbitals, 
and the dominant next nearest neighbor hoppings lead to 
quasi-one-dimensional bands, which give relatively strong 
nesting of the Fermi surface. Considering the quasi-one- 
dimcnsional as well as the doped-band-insulator nature 
of the band structure, there are several interesting pos- 
sibilities of the superconducting state depending on the 
dominating many-body interactions. 
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spin-fluctuation-mediated pairing interaction, we find the 
sign reversing s-wave gap shown in Fig|6l^b) for S ~ 0.5 
and U — 2.43eV and the c?-wave superconducting gap 
shown in Figllfc) S — 0.25 and U = 1.8eV. These gaps 
are formed by cooperation/competition between the re- 
pulsive (sign reversing) spin-singlet pairing interactions 
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